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Abstract: 

We calculate the 0(a) electroweak corrections to charged- and neutral-current deep- 
inelastic neutrino scattering off an isoscalar target. The full one-loop-corrected cross 
sections, including hard photonic corrections, are evaluated and compared to an earlier 
result which was used in the NuTeV analysis. In particular, we compare results that 
differ in input-parameter scheme, treatment of real photon radiation and factorization 
scheme. The associated shifts in the theoretical prediction for the ratio of neutral- and 
charged- current cross sections can be larger than the experimental accuracy of the NuTeV 
result. 



October 2003 



1 Introduction 



Deep-inelastic neutrino scattering has been analyzed in the NuTeV experiment [[Q with 
a rather high precision. In detail, the neutral- to charged-current cross-section ratios [|2] 

have been measured to an accuracy of about 0.2% and 0.4%, respectively. In addition, 
the quantity 

^cc(^ - A*"*) - ^gc(^ - ' 
as proposed by Paschos and Wolfenstein [Oj, has been considered. As a central result, 
the NuTeV collaboration has translated their measurement into a value for the on-shell 
weak mixing angle, sin 2 9w = 1 — M^/M^, which can, thus, be viewed as an independent 
(but rather indirect) determination of the W- to Z-boson mass ratio. The NuTeV result 
on sin 2 #w is, however, about 3a away from the result obtained from the global fit [Ej of 
the Standard Model (SM) to the electroweak precision data. 

The 3cr difference of the NuTeV result to the SM global fit to data, together with 
the fact that the evaluation of the NuTeV measurement is very involved, triggered some 
discussion about the reliability of the NuTeV result. A detailed discussion of various 
experimental as well as theoretical issues can be found in Ref. where in particular 
QCD corrections (see also Ref. [E] and references therein), charm-mass effects, and un- 
certainties from the parton distribution functions have been considered. According to 
Ref. [Uj, none of these issues seems to be a candidate for an error or an underestimate 
of uncertainties in the analysis. However, the authors of Ref. pointed out that the 
inclusion of electroweak radiative corrections, which influences the result significantly, is 
based on a single calculation [|7j only 1 and that a careful recalculation of these corrections 
is desirable. 

In this paper we present such a calculation of the 0(a) electroweak corrections to 
neutral- and charged-current deep-inelastic neutrino scattering off an isoscalar target. 
Apart from a different set of parton densities and input parameters the most important 
difference between our calculation and the result of Ref. [ Ej lies in the treatment of initial- 
state mass singularities due to collinear radiation of a photon from the incoming quark. 

The input parameters used for the numerical analysis of Ref. [Ej are not completely 
specified and, in order to fix all free parameters of the electroweak SM, need to be sup- 
plemented by fermion masses that we obtained from references quoted within Ref. [|7j. A 
comparison of numerical results should therefore be interpreted with care. In any case, 
the effect on the electroweak corrections on the NuTeV results should be further ana- 
lyzed. We will provide a Fortran code that could be used for an update of the NuTeV 
data analysis. 

The paper is organized as follows. In Section |21 we set our conventions, define the 
necessary kinematical quantities, and consider the lowest-order cross sections. Section |31 

1 Electroweak radiative corrections to neutral-current-induced neutrino deep- inelastic scattering were 
also investigated in Ref. [|Hj, but the corresponding hadronic cross sections were not evaluated numerically. 
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Figure 1: Lowest-order diagrams for v^q — > v^q and v^q — > // q' scattering 

contains an overview over the structure of the calculated electroweak corrections, a de- 
scription of the calculational framework, and explicit results on some specific one-loop 
corrections; particular attention is paid to the treatment of collinear photon radiation and 
the related fermion-mass singularities. Our numerical results are presented in Section 0] 
Section contains a short summary. 

2 Conventions and lowest-order cross sections 

We consider the neutral-current (NC) and charged-current (CC) parton processes 

NC: Vy.{pi) +q(p q ) -> v^{ki) + q{k q ), q = u, d, s, c, u, d, s, c, (2.1) 
CC: Vfiipi) + q{p q ) -> H~(ki) + q'(k q ), q = d, s, u, c, q' = u, c,d,s, (2.2) 

and the processes with all particles replaced by their antiparticles. The assignment of 
momenta is indicated in parentheses. At lowest order the former processes proceed via Z- 
boson exchange, the latter via W-boson exchange; the corresponding Feynman diagrams 
are shown in Figure ^ The diagrams and results for incoming antiquarks can be obtained 
from those of incoming quarks by crossing the external quark lines. CP symmetry implies 
that the (parton) cross sections for z/ M -quark and i/ M -antiquark scattering are equal to the 
ones of ^-antiquark and z/^-quark scattering, respectively. In the following we neglect 
the muon mass and the parton masses whenever possible. The quark masses are only kept 
as regulators for mass singularities, which appear as mass logarithms in the (photonic) 
corrections; the issue of quark-mass logarithms will be discussed below in detail. The 
muon mass is not only required as regulator for collinear final-state radiation but also for 
a proper description of forward scattering in the (loop-induced) ^v^v^ vertex, as described 
below as well. The Mandelstam variables of the partonic processes are defined as 

s = (pi+ Pq ) 2 = AE 2 , 

t = {pi-hf = -2£ 2 (l-cos#), 

u = ( Pq - ktf = -2£ 2 (l + cosfl), (2.3) 

where E and 9 denote the beam energy and scattering angle in the partonic centre-of-mass 
(CM) frame. Moreover, we introduce the variable 

y = -- = \{l-<x*0). (2.4) 
s 2 

The momentum of the incoming (anti-)quark q is related to the total nucleon momen- 
tum Ppf by the usual scaling relation 

< = xP&, (2.5) 
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which holds in the CM frame of the partonic system where the nucleoli mass is 
negligible. The variable x is the usual momentum fraction, restricted by < x < 1. 
Since deep-inelastic neutrino scattering usually is accessed in fixed-target experiments (as 
NuTeV), we take the energy E^ AB of the incoming (anti-)neutrino beam to define the 
incoming momentum p\ in the rest frame of the nucleon, called LAB in the following. 
Thus, the scattering energy E and the variable s are given by 



AE 2 = 2xM N E}l AB 



(2.6) 



up to terms of higher order in the nucleon mass. Neglecting the nucleon mass, there is also 
a simple relation between the scattering angle 9 in the partonic CM frame (or equivalently 
y) and the energy E\^ AB of the outgoing quark in the LAB frame, 



LAB 
had 



E 



LAB 



COS ( 



LAB 



(2.7) 



i^had can be identified as the energy deposited by hadrons in the detector. More generally, 
a momentum with components in the CM frame receives the following components in 
the LAB system, 



-"■LAB 



Oy(K° + (3K 3 ), K 1 , K 2 , j{(3K° + K [ 
with 7 ~ Je^ ab /(2xMj 



'iVJ 



1. 



(2- 



Note that the asymptotic expressions for 7 and f3 are fully sufficient and can be consistently 
used; deviations from the full Lorentz boost are suppressed by nucleon-mass terms. 

The lowest-order matrix elements corresponding to the tree diagrams of Figure Q can 
be easily worked out; the result is 



^NC.O 



■Mcc.0 



e2 9lqz9u^Z 

t-Ml 

2 — — 
e 9q' q w9v^fjW 

t-M& 



u(/ci)7 P ^u(p ; )J \u(k q )j p u T u(p q ) 
u(k) Y uJ-u(pi) u{k ql )^ p u_u{p q 



(2.9) 



where u>± = ~(1 ± 75) denote the chirality projectors and an obvious notation for the 
Dirac spinors u(p q ), etc., is used. The index r = ± indicates the chirality of the incoming 



quark. The coupling factors are given by 



9q'qW 



y/2. 



9lalJW 



V2& 



5 T -. 



s J? 

9hz = -fQf + r-jr 



6 T - 



(2.10) 



where Qf and J| = ±1/2 are the relative charge and the third component of the weak 
isospin of fermion /, respectively. The weak mixing angle is fixed by the mass ratio 
Mw/Mz, according to the on-shell condition sin 2 6w = s^y = 1 — c 2 ^ = 1 — .\ / 2 V ,' Mj. Note 
that the CKM matrix element for the q'q transition, V q i q , appears only as global factor 
|y g ' g | 2 in the CC cross section. Thus, in the limit of vanishing final-state quark masses, 
the factors \V q > q \ 2 add up to one after the sum over the flavours of q' is taken. 
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3 Calculation of radiative corrections 

3.1 Overview and calculational framework 

We have calculated the cross sections to the processes 1)2.1)1 and ()2.2)1 including relative 
0(a) electroweak corrections. This means our partonic cross sections a can be written as 
a sum of leading-order (er ), virtual (oVirt), and real (cr rea i) contributions: 

a — a + er virt + CT rea l = cr (l + 5), (3.1) 

where 6 denotes the relative correction. The virtual correction <j v i rt contains one-loop 
diagrams and counterterms, 

&virt = ^one-loop + °"ct- (3-2) 

The generic diagrams for these contributions, with loops and counterterms summarized 
by the blobs, are shown in Figure El Within the real correction o" rea i we can discriminate 
bremsstrahlung corrections (o"brem) and a counterterm contribution (cr co u) that accounts 
for the absorption of initial-state mass singularities into the parton densities, 

C"real = C"brem + O"coll- (3.3) 

The bremsstrahlung diagrams are shown in Figure 01 

The contributions to cr virt and cx real diverge individually and cannot be evaluated unless 
suitable regularization and renormalization prescriptions are employed. All parts of our 
calculation have been performed in two independent ways, resulting in two completely 
independent computer codes. Both loop calculations are carried out in 't Hooft-Feynman 
gauge and are based on the standard techniques for one- loop integrations as, e.g., de- 
scribed in Refs. [ El El] • Ultraviolet divergences are treated in dimensional regularization 
and eliminated using the on-shell renormalization scheme [ [HB E] i n the formulation of 
Ref. [EH|. Infrared (i.e. soft and collinear) singularities are regularized by an infinitesimal 
photon mass and small fermion masses. The artificial photon-mass dependence of the 
virtual and (soft) real corrections cancels in the sum of both contributions, according 
to Bloch and Nordsieck [ ^] . The role of the other fermion-mass singularities will be 
discussed below. 
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Figure 3: Bremsstrahlung diagrams for the processes v^q 
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v^q and v^q -> // g'. 



In detail, the two calculations of the 0(a) corrections have been performed as follows: 

1. In the first calculation, the fully differential partonic 2 — > 2, 3 particle cross sections 
have been obtained with the Mathematica packages FeynArts [ El and Form- 
Calc [ 14J. The one-loop integrals are evaluated with the LoopTools library [ 
ITlj . The combination of virtual and real corrections is done by phase-space slic- 
ing in the soft-photon region, i.e. soft photons are excluded from the phase-space 
integral of the bremsstrahlung by a small energy cutoff AE < E 1 and added as 
correction factor to the lowest-order cross section derived from the eikonal current 
(see e.g. Refs. [EI1E])- I n the whole calculation, fermion masses are kept small 
but non-zero in the matrix elements and in the phase space. 

2. The second calculation is completely independent of the first. The results for the CC 
processes f!2.2|) are obtained via crossing the results of Ref. [|T3] where the related 
Drell-Yan-like process qq' — > l~T>\ was treated in 0(a). Of course, the W-boson 
width, needed there to describe the W-boson resonance, is set to zero in our case. 
The corrections to the NC processes (|2.ip have been calculated in the same way as 
described in Ref. [ US] for the CC case. This means that the one- loop amplitudes are 
generated with FeynArts [|T3| and algebraically reduced using FeynCalc [HH] or 
own independent Mathematica routines; the scalar and tensor one-loop integrals 
are evaluated using the methods and results of Refs. [HE QUI- The bremsstrahlung 
amplitudes are derived in the Weyl-van der Waerden spinor technique, following 
the formulation of Ref. [Ej. The combination of virtual and real corrections is 
performed in three different ways: using two different versions of phase-space slicing, 
which are described in Ref. [|T3] in detail, and the dipole subtraction approach, as 
described in Ref. [El- Using these techniques, masses of the external fermions can 
be set to zero in all bremsstrahlung amplitudes consistently; the mass dependence 
relevant for mass-singular regions is restored according to factorization properties 
of amplitudes. 

We have checked that the different calculations lead to results on the relative 0(a) cor- 
rections to the NC and CC cross sections that typically agree within ~ 0.1%. For the 
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Figure 4: Diagrams for the one-loop correction to the jv^f^ vertex. 

0(a)-induced shift Asin 2 ^ in the weak mixing angle (see Section I4~T1 below) we find 
agreement within ~ 0.0003. The residual differences, which are due to the various ways 
of treating the fermion-mass effects, are thus below the NuTeV theoretical uncertainties 
related to other effects [EJE|. 

In the following, we do not present the full formulas for all contributions to the cross 
section 2 , since their derivation is by now a standard exercise and their form does not reveal 
further insight. We restrict ourselves to sketching salient features and delicate parts. 



3.2 Contributions from 2-particle final states 
tions 



tree level and virtual correc- 



The differential cross section including all 2 
given by 

da 1 1 E(l-^5l 2 + 2 M(-M5r-MI} 



dcos# 2 327rs 



2 particle contributions at 0(a 3 ) is 

(3.4) 



where M.q and A4\ are the amplitudes at tree and one-loop level, respectively. Note that 
rotational invariance has been used to integrate over the azimuthal angle </>, leading to 
a global factor 2ir, and the explicit prefactor 1/2 accounts for the spin average of the 
initial-state quark. 

Among the virtual corrections, there are two contributions to the NC processes that 
become numerically delicate in the limit of small momentum transfer (t — > 0): the ^v^v^ 
vertex correction and the 7Z mixing self-energy, indicated in the last two generic diagrams 
in the first line of Figure El The corresponding contributions to the relative correction S 
read 



5 T - 

7»W 



4s w c w Q g t - Mf 



9 T qqZ 



(*) + 



As 



-8Z 



2Q q 

9qqZ 



+ -5Z lZ + 



t-Ml 



2t 



-5Z 



(3.5) 
(3.6) 



where is the transversal part of the unrenormalized 7Z mixing self-energy and bZ^z 
and bZz^ are the renormalization constants for the mixing of the 7 and Z fields. The 
precise definition and explicit expressions for these quantities in 't Hooft-Feynman gauge 
can be found in Ref. [CHj- The one- loop diagrams contributing to the 7^^ vertex are 
depicted in Figure where <fi represents the would-be Goldstone field corresponding to the 
W boson. The renormalized vertex form factor F™ p (t), which is given by the expression 

2 For the CC processes all contributions can be obtained from the formulas of Ref. [El y i a crossing. 
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in square brackets in Eq. (|3.5|) . explicitly reads 

a 



F Ten - 



(*) 



16ns 2 



2(m 2 + 2M 2 ) 

^ W - 4[5 (0, M w , M w ) - S (t, M w , M w )] 



2m* + 2m 2 M w - 4M W + m 2 / - 6M W / 



£M W 



x [B (t, m M , m M ) - B„(i, M w , M w )] 

2(m^ - 3m 2 M w + 2M W + mfo - 2m 2 M w t + 4M w t) 



+ 



*M W 



x [C (0, i, 0, m M , M w , M w ) + C Q (0, t, 0, M w , m M , m M )] 

+ 4tCo(0,t,0,M w ,m M ,m M )|. (3.7) 

The definitions and general results for the scalar 2- and 3-point integrals Bq and Co can 
be found in Ref. [HUj. Note that the full muon mass dependence is kept in this formula. 
For large \t\, i.e. for \t\ 3> m 2 , the muon mass can be neglected; in this limit the form 
factor reads 



rpren 

7^M P M 



(*) 



a 



+ 



167TS 2 

2(2M W + 3/) 
t 



4 - 4[5 (0, M w , M w ) - B {t, M w , M w )] 



;5 (t,0,0)-Bo(t,M w ,M w )] 



4 /If 2 I'M 2 -I- 2/1 

+ wl , w+ J [C (0, /, 0, 0, M w , M w ) + C (0, /, 0, M w , 0, 0)] 



+ 4/C (0,/,0,M w ,0,0) 



(3i 



which reproduces the full expression for up to terms of C(m 2 //) and C(m 2 /M w ). 

For small |/|, the general expression ()3.7|) becomes numerically unstable, because the 
whole expression in this limit is of 0(t) while individual terms involve constants and even 
\jt poles. A stable result for this limit can be obtained by performing an asymptotic 
expansion for large Mw; this leads to 



/? ren 



(*) 



a 



6m 2 [B (0, m„, m M ) - B (t, m M , m M )] 



Mw-oo 367rM^s 2 v 

*[2 - 3,00(0, M w , M w ) + 3B (t, m M )] k 



(3.9) 



which is valid up to terms of 0(mJ/M w ), C(m 2 //M w ), and C(/ 2 /M w ). From this ex- 
pression one can easily obtain the limit £ — > for the contribution of the ^v^v^ vertex to 
the relative correction, 



7^M P M 



Qo 



a 



t=0 



(3.10) 



In the limit t — > 0, the relative correction induced by 7Z mixing tends to the value 



t=o 



2Q q 

9qqZ 



E^ Z (0) -S^(M| 



+ (0) 



M 2 



(3.11) 



where T,^ z '(t) is the derivative of E^ z (t) w.r.t. the variable t. Since the derivative S^ Z/ (0) 
numerically results from [E^ z (t)— Eij Z (0)]/t for t — > 0, the actual evaluation of deserves 
some care for small \t\. This is, however, less tricky than for 8^ v p , since only one power 

of £ effectively cancels for £ — > 0. A peculiarity in E^ Z/ (0) is that it involves logarithms 
ln(m//Mz) of all charged fermions /. These contributions result from the jZ mixing in 
the non-perturbative domain of small momentum transfer. Thus, the light-quark masses 
entering here should be tuned in such a way that the hadronic jZ mixing, which can 
be obtained via dispersion relations is reproduced. This highly non-trivial task is, 

however, beyond the scope of this paper. 

In the actual numerical evaluation, we proceed in the two following ways (correspond- 
ing to the two independent calculations described under the same labels in Section l3~T|) : 



1. In the first calculation we substitute 5Z„ 



and 51 



expression if |£| < 5 x 1CT 3 GeV 2 . 



for the exact i-dependent 



2. In the second calculation the asymptotic formulas ()3.8|) or ()3.9|) are used if |£| is larger 
or smaller than lGeV 2 , respectively. At the matching point, these two formulas 
agree within about 4 digits. The ^Z mixing part is evaluated with the full expression 
(|3.5|) everywhere. 

3.3 Contributions from 3-particle final states — bremsstrahlung corrections 

For the real photonic corrections, the processes (J2.1|) and (|2.2j) have to be considered 
with an additional photon (with momentum k) in the final state. The contribution Ubrem 
of these radiative processes to the parton cross section are given by 



1 1 

<Xbrcm - 2 



J dr 7 £ \M,\\ (3.12) 



spins 

where M. 1 is the corresponding amplitude and the phase-space integral is defined by 
d 3 ki f d 3 k q f d 3 k 



ir -- = J wfki J J M ^ + p ' - h ~ k ° ~ ky (X13) 

The phase-space integral ()3.12|) diverges in the soft (k — > 0) and collinear (p q k, k q k — > 0, 
and k\k —>■ in the CC case) regions logarithmically if the photon and fermion masses 
are set to zero. For the treatment of the soft and collinear singularities we have applied 
the phase-space slicing method as well as a subtraction approach, as mentioned above. 

3.4 Hadron cross section and redefinition of parton densities 

The neutrino-nucleus scattering cross section <t u n is obtained from the parton cross 
sections o by convolution with the corresponding parton distribution functions q(x), 

a uN (K AB ) = E [ dxq iso {x)a{p hPg ), (3.14) 

where x is the momentum fraction carried by the parton q and M is the factorization 
scale. The squared CM energy s of the partonic process is related to the neutrino energy 
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^lab j n frame according to Eq. ()2.6|) . In the sum J2 q the variable g runs over the 

incoming quarks and antiquarks indicated in Eqs. (|2.1|) and (|2.2j) . The subscript "iso" 
indicates that we average over the parton densities for an (anti-) quark q in the proton 
and neutron, q(x) = f q (x) and f q (x), 



(3.15) 



Owing to isospin invariance (f^(x) = f^(x), etc.), this is equivalent to an isospin average 
of the parton densities of the proton, u iso (x) = d iso (a;) = |(u(a;) + d(x)), etc. 

The 0(a)-corrected parton cross section o contains mass singularities of the form 
a\n(m q ), which are due to collinear photon radiation off the initial-state quarks. In com- 
plete analogy to the MS factorization scheme for next-to-leading order QCD corrections, 
we absorb these collinear singularities into the quark distributions. This is achieved by 
replacing q(x) in (j3.14J) according to 



q(x) q(x,M 2 ) 



P, 



ff\ 



P ff (z)(2\n(l-z) + l 



(3.16) 



where M is the factorization scale (see Ref. [E])- This replacement defines the same finite 
parts in the 0(a) correction as the usual MS factorization in D-dimensional regularization 
for exactly massless partons, where the ln(m g ) terms appear as l/(D — 4) poles. In (|3.1fij) 
we have regularized the soft-photon pole by using the [...]+ prescription. This procedure 
is fully equivalent to the application of a soft-photon cutoff AE (see Ref. [1201) where 



q(x) 



q(x,M 2 



1 _ ?L q4\ _ ln(2A J E/v / S) - ln 2 (2A J E/v / S) 



1-2AE/VJ dz 

Z 



M 2 



ml (1 — z) 2 1 



1 . (3.17) 



The convolution of the correction terms to the parton densities [r.h.s. of Eqs. (|3.1fij) or 
()3.17j) ] with the lowest-order cross section <7o defines the "collinear counterterm" contri- 
bution cr co n in the real photonic correction <7 re ai, as defined in Eq. ()3.3|) . 

The absorption of the collinear singularities of 0(a) into quark distributions, as a 
matter of fact, requires also the inclusion of the corresponding 0(a) corrections into 
the DGLAP evolution of these distributions and into their fit to experimental data. At 
present, this full incorporation of 0(a) effects in the determination of the quark distri- 
butions has not yet been performed. However, an approximate inclusion of the 0(a) 
corrections to the DGLAP evolution shows [|^] that the impact of these corrections on 
the quark distributions in the MS factorization scheme is smaller than remaining QCD 
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uncertainties; this is also supported by a recent analysis of the MRST collaboration [ |2*2~] 
who took into account the 0{a) effects to the DGLAP equations. 3 

3.5 Structure of divergences in the radiative corrections 

While ultraviolet divergences are absorbed by renormalization and soft singularities 
always compensate between virtual and (soft) real corrections, the issue of collinear sin- 
gularities is more delicate. These singularities originate from collinear photon emission 
off the external fermions and lead to mass logarithms alnrrif. As described in the pre- 
vious section, the a In m q terms from initial-state radiation are absorbed into the parton 
distribution functions. 

However, there are also a In m q > and (in the CC case) a In m M terms from final-state 
radiation. According to the Kinoshita-Lee-Nauenberg (KLN) theorem these terms 
drop out if the final state is treated sufficiently inclusive, i.e. if the cones for quasi-collinear 
photon emission around the charged outgoing fermions are fully integrated over. This con- 
dition is, in general, not fulfilled if phase-space cuts at the parton level are applied. For 
instance, restricting the outgoing hadronic energy, E^^ min < E^^P < E^^ max , spoils the 
inclusiveness w.r.t. photon emission collinear to the outgoing quark, so that corrections 
proportional to a \nm q i survive. On the other hand, a cut on the total hadronic+photonic 
energy, £had+ P hot,min < ^had+phot < £had+ P hot,max, which is equivalent to cutting the out- 
going lepton energy, destroys the inclusiveness w.r.t. collinear photon emission of muons 
in the CC processes, leading to a In corrections. Both types of mass singularities are 
avoided if photon recombination in the final state is taken into account before applying 
selection cuts. 

For instance, a simple recombination scheme is provided by the following algorithm: 

1. For an event with a photon in the final state, the smallest angle 9 7 f between the 
outgoing photon and the outgoing fermions (q for NC; \i or q' for CC) is determined 
in the LAB frame. 

2. If 6 1 f < 6 cut, the photon is recombined with the nearest fermion /. This means that 
the momenta of the photon and of the fermion / are added and assigned to the new 
"quasi-fermion" /; then the photon is discarded from the event. 

3. Finally, we impose the hadronic energy cut ^had B reC omb > 10 GeV, where "recomb" 
indicates that the momentum after possible recombination is used. 

This procedure is "collinear safe" with respect to initial-state and final-state collinear 
singularities, i.e. photons that are nearly collinear to a charged fermion are treated inclu- 
sively. Note that photons exactly collinear to the initial-state quark receive zero momen- 
tum in the LAB frame, which effectively acts as recombination with the initial state. All 
final-state muon and quark-mass logarithms cancel in the corresponding integrated cross 
section. The collinear safety of observables can be ideally checked by comparing results 

3 Note that the inclusion of 0(a) effects in the DGLAP evolution also leads to a photon distribu- 
tion function in the nucleon. The impact of the resulting v^j reactions, which are connected to the 
bremsstrahlung corrections via crossing, is, however, expected to be very small owing to smallness of the 
induced photon distribution. 
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obtained with phase-space slicing to the ones obtained with a subtraction technique; we 
have performed this comparison and find agreement within the integration errors. In the 
numerical evaluation we take 8 cnt = 5°. 

Although all of the above cuts are idealizations, the cut on the sum -Ehad+phot °f 
hadronic and photonic energies seems to be most appropriate for a fixed-target experiment 
where all but the energy of neutrinos and muons is deposited in a calorimeter. 

4 Discussion of results 

4.1 The quantities 8R U and A sin 2 W 

Although the electroweak 0(a) corrections to the NC and CC cross sections represent 
the central results of this paper, it is useful to translate them into a correction to the 
cross-section ratio R v = Onc/^cc introduced in Eq. and to make contact to its 

interpretation in terms of an effective shift in the on-shell weak mixing angle. 

Writing 5a for the correction to the cross section, the (first-order) correction 8R V to 
the ratio R v reads 

Sit = R» ( ^ - R" (SR^ C + 5R» CC ) . (4.1) 



a NC a CC 



Following Ref. [|7], the quantities 

ff NC (T CC 

have been defined for later convenience. 

In order to define an effective shift Asin 2 6\y, we quote the approximate relation 

1 20 

which results from the lowest-order neutrino deep-inelastic scattering cross section ratio 
for isoscalar targets upon neglecting the momentum transfers in the W- and Z-boson 
propagators. The correction 5R U can be interpreted as a correction 5 sin 2 #w in the variable 
sin 2 W , 

5R V ~ S Q - sin 2 6 W + P sin 4 W ) = (-1 + — sin 2 9^) 6 sin 2 9 W . (4.4) 

As done in Ref. [ |7] , we define the shift A sin 2 6*w as the negative of the correction S sin 2 #w- 
Combining the above relations, we are, thus, lead to 



A sin V w = -d sin W = 40 2fl — — — 



k - sin 2 6>w + if sin 4 8 W ( Sa^ c 5a l c 

40 2 n \ -v -v 

27 sin V w \ a NC a cc j 

|-sin 2 fl w + §sin 4 fl w , v , , 

40 -2n ( dK NC + dK CC> • I 4 - 5 ) 

k= sm t/w 



27 
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4.2 Input parameters 

If not stated otherwise, we used the following set of input parameters [ Olj for the 
numerical evaluation, 





= 1.16639 x Kr 5 GeV" 2 , 


a(0) 


= 1/137.03599976, 


a(Afz) 


= 1/128.930, 


M N 


= m p = 0.938271998 GeV, 










M w 


= 80.423 GeV, 


M z 


= 91.1876 GeV, 


M H 


= 115 GeV, 


m e 


= 0.510998902 MeV, 


rrift 


= 105.658357 MeV, 


m T 


= 1.77699 GeV, 


m u 


= 66 MeV, 


m c 


= 1.2 GeV, 


m t 


= 174.3 GeV, 


m d 


= 66 MeV, 


m s 


= 150 MeV, 


m h 


= 4.3 GeV. 



(4.6) 

The masses of the light quarks are adjusted such as to reproduce the hadronic contribution 
to the photonic vacuum polarization of Ref . [ • They are relevant only for the evaluation 
of the charge renormalization constant in the a(0)-scheme and for the small correction 
induced by the jZ mixing at zero- momentum transfer [see Eq. (J3.11|) ]. For the calculation 
of the u^N cross sections we have set the energy of the incoming neutrino to E^ AB = 
80 GeV and adopted the CTEQ4L parton distribution functions [:.26]- The factorization 



scale M is derived from the leptonic momentum transfer, i.e. it is set to M = y—(pi — ki) 2 
(with the momenta after possible photon recombination). 

For the explicit evaluation and the following numerical discussion we distinguish three 
different input-parameter schemes 

• a(0)-scheme: The fine-structure constant a(0) and all particle masses define the 
complete input (i.e. a(Mz) and G M are not used). In this scheme, the relative 
corrections to the cross sections sensitively depend on the light quark masses via 
a\nm q terms that enter the charge renormalization. 

• a(Mz) -scheme: The effective electromagnetic coupling a(Mz) and all particle 
masses define the complete input (i.e. a(0) and G^ are not used). Tree-level 
couplings are derived from a(Mz), and the relative corrections are related by 
<5a(Af z )-scheme = ^(o)-schcme - 2Aa(M z ), where Aa{Q) accounts for the running 
of the electromagnetic coupling from scale Q = to Q = M z (induced by light 
fermions) and cancels the corresponding a\nm q terms in £ a (o)-scheme- 

• G ^-scheme: The Fermi constant G^ and all particle masses define the complete 
input (i.e. a(0) and a(M z ) are not used). Tree-level couplings are derived from the 
effective coupling «g m = y/^G^My^s^/iT, and the relative corrections are related by 
<5g m -scheme = <5 a (o)-scheme — 2 Ar , where Ar contains the radiative corrections to muon 
decay and is defined [EZj by 

M 2 ( 1 _Mk) = _12L \ a 7) 

w \ Mi) ^2G,l-Ar(a,M w ,M z ,M K ,m f Y { ' } 

In the following we consistently use Ar at the one-loop level, as, for instance, ex- 
plicitly given in Ref. [E3]. Since Aa(Mz) is contained in Ar, there is no large effect 
induced by the running of the electromagnetic coupling in the G M -scheme either. 

More details about the various schemes can, e.g., be found in Refs. [ FTTH ITU ITB] . 



12 



Hadronic energy cut: E^Qf > 10 GeV 



IPS 


Rq 


i?ATn 


5R V CC 


A sin 2 #w 


a(0) 

Qt(Mz) 


0.31766(2) 
0.31766(2) 
0.31766(2) 


0.0582(1) 
-0.0639(1) 
0.0003(1) 


-0.0758(4) 
0.0452(4) 
-0.0185(4) 


-0.0082(2) 
-0.0088(2) 
-0.0085(2) 


Hadronic+photonic 


energy cut: 


^ha A d B + phot> 10GeV 


IPS 






5R? CC 


A sin 2 # w 


a(0) 
a(Af z ) 


0.31766(2) 
0.31766(2) 
0.31766(2) 


0.0589(1) 
-0.0632(1) 
0.0011(1) 


-0.0842(4) 
0.0363(4) 
-0.0272(4) 


-0.0118(2) 
-0.0126(2) 
-0.0122(2) 



Hadronic energy cut after 7 recombination: ^hadrecomb > 10 GeV 



IPS 


Rq 


5R^ C 


5R V CC 


A sin 2 6'w 


a(0) 


0.31766(2) 


0.05861(1) 


-0.07702(2) 


-0.00863(1) 


a(M z ) 


0.31766(2) 


-0.06349(1) 


0.04392(2) 


-0.00917(1) 




0.31766(2) 


0.00077(2) 


-0.01980(2) 


-0.00892(1) 



Table 1: Results on the ratio R v in lowest order (Rq), corrections from NC and CC 
cross sections (SR^ C and 5Rq C ), and shift A sin 2 6'w, for various input-parameter schemes 
(IPS) and using different treatments of real photons in the final state. The number in 
parentheses indicates the statistical integration error in the last digit. 



4.3 Numerical results 

Using the input defined in the previous section, we have evaluated the corrections SR^ C 
and SRq C , as well as the corresponding shift A sin 2 6'w, in the various input-parameter 
schemes and for different treatments of photons in the final state. The corresponding 
numerical results are collected in Table Q 

We observe a large scheme dependence of the corrections SR^ C and SRq C owing to the 
shifts of — 2Aa(Mz) ~ —12% and — 2Ar ~ —6% in the relative corrections 5 Q (M z )-scheme 
and ^G^-scheme, respectively, relative to the a(0)-scheme. Since these shifts are identical 
in the NC and CC cross sections, they drop out in A sin 2 6'w- The remaining very small 
scheme dependence in A sin 2 6'w is due to the different choice of the coupling a in the 
schemes, i.e. of formal two-loop order. 

Moreover, Table reveals that different ways of treating real photons in the final state 
selection lead to sizeable differences in the correction 5Rq C to the CC cross section while 
5R^ C hardly changes. The differences can be explained by inspecting final-state radiation 
w.r.t. their inclusiveness. The variant -^had+phot > 10 GeV, where the energy cut is applied 
to the sum of hadronic and photonic energy, is the only choice in which the collinear cone 
for radiation off outgoing muons is not treated inclusively, leading to logarithmically 
enhanced corrections proportional to alnm^. Numerically this effect changes SRq C by 
~ 0.008, and thus increases the shift A sin 2 6'w by ~ 0.004. On the other hand, the variant 
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energy cut is applied to the bare outgoing quark energy, is the 
only case in which final-state quark logarithms alnm q (in the NC case) or a\nm q i (in 
the CC case) are present. Numerically these effects turn out to be of the order ~ 0.001 
in an d <^cc> anc ^ thus do not influence A sin 2 #w very much. 

4.4 Comparison with existing results 

Finally, we compare our results with the ones of the older calculation of Ref. [Ej as 
far as possible. To this end, we adopt the (obsolete) input parameters of Ref. [ E| which 
are, however, unfortunately not completely specified there. 4 The following parameters are 
explicitly given: 

a(0) = 1/137.036, = (37.281 GeV) 2 , sin 2 6 W = 0.227, 

v2G M 

M z = 93.8 GeV, M H = 100 GeV, m t = 180 GeV. (4.8) 

Neither the W-boson mass nor the masses of fermions other than the top quark are given 
in Ref. [ E| . We supplement the above input parameters by the following fermion masses, 

M N = m p = 0.938271998 GeV, 

m e = 0.5110034 MeV, = 0.10565943 GeV, m T = 1.7842 GeV, 

m u = 80 MeV, m c = 1.5 GeV, 

m A = 80 MeV, m s = 0.30 GeV, m h = 4.5 GeV. (4.9) 



In detail, the lepton masses are taken from Ref. [ |28] and the quark masses from Ref. 
Ref. [Ej indirectly refers to these papers. 

The authors of Ref. [Ej advocate the on-mass-shell renormalization scheme with the 
fine-structure constant a(0), the Fermi constant G^, the Z- and Higgs-boson masses Mz 
and Mh, and all fermion masses as physical input parameters. In contrast to the conven- 
tional on-shell scheme the W-boson mass Mw is not an input parameter but is instead 
determined by (iterative) solution of the implicit relation ([4.7)1 . Using the above input 
parameters, except for sin 2 9 W , and solving Eq. ()4.7|) for M w we find for the W-boson 
mass Mw = 83.807 GeV, which leads to s 2 ^ = sin 2 = 0.2017. This constitutes roughly 
a 10% discrepancy to the value sin 2 = 0.227 given in Ref. [Ej; it is not clear where 
the latter is actually used in Ref. [|7j. Having specified Mw, we further assume that the 
authors of Ref. [|7j proceed as in our G M scheme, i.e. by taking G^ and all particle masses 
as input. 

Concerning the treatment of initial-state quark-mass singularities, there is a conceptual 
difference between the calculation of Ref. [Ej and ours. The authors of Ref. [Ej do not 
apply MS factorization, as we do, but instead set the masses m q of the incoming quarks to 
m q = xMjq without any additional subtractions. Moreover, the (rather old set of) parton 
distribution functions used in Ref. [El are not available to us; the corrections, however, 
should depend on the parton densities only very weakly. 

In Table |2] we compare our results with the ones quoted in Ref. [ E] ■ Since the input is 
not completely clear, we use different setups. Specifically, we apply the usual G^-scheme 

4 The full set of input parameters used in Ref. [Cj could also not be provided by the authors on request. 



14 



Hadronic energy cut: -Ehad > 10 GeV 



IPS 


FS 


Rq 


5R^ C 


5R V CC 


A sin 2 #w 


result of Ref. [U\ 


BD 




-0.0021 


-0.0223 


-0.0114 


Gfi 


MS 


0.31455(1) 


0.0010(1) 


-0.0202(4) 


-0.0090(2) 


G M (M W derived) 


MS 


0.33113(2) 


-0.0018(1) 


-0.0186(4) 


-0.0095(2) 




BD 


0.31455(1) 


-0.0026(1) 


-0.0184(4) 


-0.0098(2) 


G (i (Mw derived) 


BD 


0.33113(2) 


-0.0050(1) 


-0.0169(4) 


-0.0103(2) 



Hadronic+photonic energy cut: E\^ +phot > 10 GeV 



IPS 


FS 


Rq 


6R^ C 


5Rq C 


A sin 2 #w 


G/j, 


MS 


0.31455(1) 


0.0018(1) 


-0.0294(4) 


-0.0130(2) 


G^(Mw derived) 


MS 


0.33113(2) 


-0.0010(1) 


-0.0271(4) 


-0.0132(2) 


G> 


BD 


0.31455(1) 


-0.0018(1) 


-0.0276(4) 


-0.0138(2) 


G M (Mw derived) 


BD 


0.33113(2) 


-0.0043(1) 


-0.0253(4) 


-0.0139(2) 



Hadronic energy cut after 7 recombination: recomb > 10 GeV 



IPS 


FS 


R o 


5R^ C 


5R U CC 


A sin 2 6> w 


G/j, 


MS 


0.31455(1) 


0.00138(2) 


-0.02149(2) 


-0.00943(1) 


G^(M W derived) 


MS 


0.33113(2) 


-0.00139(2) 


-0.01970(2) 


-0.00988(1) 


Gfi 


BD 


0.31455(1) 


-0.00214(2) 


-0.01969(2) 


-0.01023(1) 


Gf,(M w derived) 


BD 


0.33113(2) 


-0.00465(2) 


-0.01805(2) 


-0.01064(1) 



Table 2: Comparison with results of Ref. [Ej. 



(denoted "G M "), where M w = 82.469 GeV is calculated from sin 2 9 W = 0.227 and M z , as 
well as a modified G M scheme [denoted "G M (Mw derived)"] that differs from the former 
in the value of Mw = 83.807, which is calculated from and the particle masses via 
Ar at one loop. Moreover, we apply two different factorization schemes, the MS scheme 
and the one used in Ref. [Cj (quoted as "BD"). Of course, this comparison of results 
can neither prove nor disprove the correctness of the results of Ref. [|7j. However, the 
first part of the table (E^!^ > 10 GeV) suggests significant differences in the correction 
A sin 2 #w I n any case, the variations in the corrections that are due to the different 
factorization schemes (MS versus BD) and due to different ways of including the final- 
state photon in the hadronic energy in the final state can be as large as the accuracy 
in the NuTeV experiment, which is about 0.0016 in sin 2 6'w (if statistical and systematic 
errors are combined quadratically) . 

5 Summary 

Results from a new calculation of electroweak 0(a) corrections to neutral- and 
charged- current neutrino deep-inelastic scattering have been presented. The salient fea- 
tures of the calculation and some delicate points of the corrections, such as the issue of 
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collinear fermion-mass singularities, have been discussed in detail. Unfortunately a tuned 
comparison to older results, which were used in the NuTeV data analysis, is not possible. 
However, a comparison based on an assumption for missing input parameters leaves room 
for speculating on possible significant differences in the electroweak radiative corrections. 
Therefore, an update of the NuTeV analysis seems to be desirable. We provide a Fortran 
code for the electroweak radiative corrections that could be used in this task. Specifically, 
our investigation of the factorization-scheme dependence for initial-state radiation and of 
different ways to treat photons in the final state reveals that these effects can be as large 
as the 3a difference between the NuTeV measurement and the Standard Model prediction 
in the on- shell weak mixing angle. 

A careful re-assessment of all theoretical uncertainties involved in the NuTeV analysis 
should be worked out, based on the new studies of both QCD [ El Ej and electroweak 
corrections. The NuTeV collaboration estimated the theoretical uncertainty due to miss- 
ing higher-order effects to 0.00005 and 0.00011 in 5R U and Asin 2 ^, respectively. The 
results on electroweak corrections presented in this paper indicate that these numbers are 
too optimistic. 
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